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Length  Distribution  of  Chords 
Through  s  Rsctsngulsr  Volums 


SlnflL«-EygD.t-Uflgatg 

This  report  presents  differential  and  integral  chard  length 
distributions  for  use  in  predicting  single  event  upsets  (SEUs)  in 
spacecraft  electronics. 

An  "upset"  is  a  change  of  state,  such  as  the  state  denoting  a 
bit  in  computer  memory.  A  "single  event"  is  an  occurrence  which 
produces  sufficient  ionization,  in  the  speck  of  silicon  involved, 
to  produce  an  upset.  Examples  of  an  event  are  alpha  emission  by 
a  radioactive  impurity  in  the  device,  the  highly-ionized  track  of 
a  heavy  cosmic  ray,  and  a  nuclear  reaction  produced  in  the  device 
by  a  high  energy  proton.  In  the  case  of  a  penetrating  particle, 
an  SEU  will  be  produced  if  (a)  the  rate  of  ionization  along  the 
path  times  <b)  the  path  length  in  the  sensitive  region  of  the 
device  exceeds  (c)  the  critical  charge  for  upset. 

Much  of  the  present  material  was  written  as  part  of  a  larger 
report.  That  report  has  been  dropped,  but  the  new  results  herein 
have  been  requested  by  a  number  of  persons  and  used  in  at  least 
four  publ i cations. 

a_v ar  gua-£haat  * 

The  distribution  of  path  lengths  of  isotropic  particles  in  a 
given  volume  is  both  a  mathematical  and  a  physical  problem.  An 
atomic  projectile  loses  energy  by  making  many  minor  collisions 
with  the  atoms  of  the  environment.  The  direction  of  travel  is 
slightly  altered  in  each  collision.  Major  collisions,  with  large 
changes  in  direction,  are  rare.  The  "pathlength" ,  or  physical 
distance  travelled,  will  exceed  the  mathematical  chord  length. 
In  this  report,  we  treat  only  the  mathematical  problem  -  the 
distribution  in  length  of  chords,  straight  lines  between  entry 
and  exit  paints,  for  penetrating  undeviating  particles  in  a 
rectangular  volume. 

Source  and. Basic  Equation 

The  chord  length  distribution  for  an  isotropic  particle  flux 
is  given,  in  differential  form,  in  an  article  by  Pickel  and 
Blandford,®  which  contains  equations  developed  by  Petroff.  <A 
Rockwell  International  internal  report  was  actually  used  here.) 
Dr.  Philip  Shapiro  of  NRL  has  confirmed  the  equations  and  then 
employed  them  in  rr_!;.»erical  integrations.*8'  The  analytic  integra- 
ManiMcript  approved  April  23,  1984.  2 


tion,  not  -found  in  the  above  source,  is  performed  below  after 
considerable  simplification  of  the  differential  form. 


We  consider  a  rectangular  parallelepiped  of  dimensions  L,  W, 
and  H.  Some  lines  will  merely  nick  a  corner  and  have  short  paths 
in  the  volume}  many  particles  (those  which  enter  a  surface  nearly 
normally  and  exit  at  the  opposite  surface)  have  a  path  a  little 
more  than  L,  W,  or  Hi  the  maximum  chord  length  is  R,  the  diagonal 
of  the  volume. 

The  normalized  differential  distribution  function  for  chord 
length  S  is  called  F  and  integrates  to  unity.  It  is  convenient 
to  replace  the  quantity  g  of  Petroff  by  M/12,  where  M  *  N  +  Mo, 
and  Mo  consists  of  those  terms  whose  contributions  are  zero.  The 
function  F  is  given  by 

3 if  CLW  +  LH  +  WHD  F(S,L,W,H>  -  2  N(S,X,V,Z),  <1> 

* 

where  the  sum  is  over  the  six  permutations  of  (L,W,H)  as  (X,Y,Z). 

As  the  equations  of  Petroff  are  quite  complex,  symbols  are 
used  here  for  certain  functions  of  the  dimensions: 

K2  .  xa  +  Y3,  R*  -  La  +  W*  +  Ha  «  X*  +  Ya  +  Z2, 

Ta  -  X2  +  Z2,  and  V  -  12  XYZ2. 

Similarly,  two  symbols  show  relationships  with  chord  length: 

P2  -  S2  -  Z2  and  Q2  -  S2  -  X2  -  Z2. 

Initial  JHffgiQtlal  JEacmu.  M 

The  expressions  far  M  and  N  have  three  forms,  depending  upon 
the  value  of  S  relative  to  the  dimensions.  First,  M  is  given  as 
a  sum  of  quantities  which  will  be  examined  individually.  In  the 
first  domain,  where  0  <  S  <  Z,  it  is  given  by 

Mi  -  N»  ♦  Moa  +  Mob.  (2) 

In  the  second  domain,  where  Z  <  S  <  T,  Petroff 's  formula  is  equal 
to 

M*  m  Na  +  Moa  ^  Mob  ^  Moc«  (3) 

Finally,  in  the  third  domain,  where  T  <  S  <  R,  the  expression  is 

M3  “  N3  +  Moa  +  Moc  (4) 

The  sum  in  Eq.  (1)  will  often  involve  different  regions  for 
different  permutati ons.  For  example,  if  L  >  S  >  H,  some  of  the 
permutations  (with  Z  ■  L)  will  involve  the  first  domain  while 
others  (with  Z  ■  H)  will  involve  the  second  or  third  domain. 

The  terms  combined  into  Moa  are  common  to  all  three  domains 
and  therefore  permutations  can  be  combined  for  all  values  of  S. 
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These  terms  are 


M0a-  8XaCYa/TaR  -  2/R  +  R/Ka3  -  <8/R> CX*Ya/Ta  +  XaZa/Ka  -X*D.  (5) 
Summing  over  all  six  permutations,  one  -finds  that 

1  Mo«  -  <8/R) C2Ra  -  2Ra3  -  0,  (Sa) 

and  Moa  may  be  omitted  (as  it  was)  -from  Eq.  (1). 

Now  let  us  consider  the  quantity 

Mo*  -  8(X2  -  Za)/T.  (6) 

The  term  appears  in  the  -first  and  second  expressions  -for  M,  or 
•for  0  <  S  <  T.  Permutations  with  Y  fixed  but  X  and  Z  exchanged 
will  have  identical  values  of  T  and  thus  the  same  range.  These 
terms  cancel  in  pairs,  so 

2  Mo*  -  0.  (6a> 

The  third  quantity  to  be  eliminated, 

Moc  -  (X*  -  Ya) (P/KS) (8  +  4Za/S2) ,  (7) 

exists  for  Z  <  S  <  R.  One  may  sum  permutations  if  Z  is  fixed. 
Again,  paired  cases  cancel  and 

I  Moc  -  0.  <7a) 

The  value  of  F  is,  thus,  obtained  from  the  sum  of  the  N’s,  as 
shown  in  Eq.  (1). 

Net  Pi  tifficmt  f grata N 

In  the  first  domain,  where  0  <  S  <  Z,  what  remains  is 

Nx  -  BiS  +  4Z(2Y=/K=),  (S> 

where  B*  ■  -(3XY/RT)3*.  It  is  easily  seen  that,  when  summed  over 
paired  permutations  of  X  and  Y,  this  is  equivalent  to 

Nik  -  BiS  +  4Z.  (8E) 

In  the  second  domain,  Z  <  S  <  T,  one  finds  that 

Na  -  BaS  +  Ba/S*  -  X (P/S) <8  +  4Za/S=>,  (?) 

where  Bs  -  (3Y/K)2  +  Bx  and  B*  -  V  tan"1 (Y/X)  -  (YZ2/K) 2. 

For  the  third  domain,  T  <  S  <  R,  the  value  is 

N*  -  B*S  +  Bs/8*  +  YZ2(Q/S>  CB/T2  +  4/S23  -  (V/S*)cos  -MX/P),  (10) 

where  3S  -  - (3XZ/KR) '  and  Ba  •  X2Z=(Z2/K2  -  3)  +  V  tan'MY/X). 


As  X  and  Z  may  be  exchanged,  this  may  bs  simplified  a  littls  to 


Na«  «  B*S  +  B»/S*  +  Y  <Q/S)  C4  +  2T*/8*3  -  <V/S*>  cos-MX/P> .  (10E) 

Short  Chord*  aori  P«iki 

Six  distributions  of  chord  lsngth  are  shown  in  Figs.  1-3. 
Ths  dimsnsions  are  in  micrometers,  as  befits  microelectronics. 
In  all  cases,  the  maximum  chord  length,  S  ■  R,  is  an  integer  or 
half -integer  from  30  to  34.  The  significant  quantities,  however, 
are  the  ratios  of  dimensions.  Defining  H  as  the  smallest,  R/H  is 
2.123,  3,  5.5,  7.5,  12.6,  and  17  for  these  cases.  The  ratio  L/W 
is  1  to  1.5. 

The  dominant  features  are  the  peaks  at  S  equal  to  H,  W,  or  L. 
Chords  of  >0.9  R  are  too  rare  to  be  evident  on  the  graphs.  It 
is  seen  that  F  is  linear  for  S  <  H.  In  this  range,  Nae  applies 
to  all  six  perturbations,  and  the  sum  is 

ZN.<h  -  8CL  +  W  +  HI  -  9S.  (11) 

The  discontinuities  can  be  evaluated  by  subtracting  NaK  from 
Na  for  the  two  permutations  involved  at  a  given  jump.  The  change 
at  S  -  Z  is  given  by 

CLW  +  LH  +  WH3  6F  -  2XY/Z.  (12) 

This  equation  shows  that  the  biggest  peak  comes  at  the  smallest 
dimension,  where  S  *  Z  *  H,  as  clearly  seen  in  the  figures.  Two 
peaks  coincide  in  a  square  device. 

Iptef  LnUU.  IntMLfllJi 

Integration  of  N1B  is  trivial!  integration  of  Na  is  quite 
easy!  integration  of  N3k  involves  a  term  not  readily  found  in  an 
integral  table.  The  results  follow: 

It  -  /n»*  dS  -  B*Sz/2  +  4ZS  (13) 

la  -  Ba85*/2  -  B*/2SZ  +  XP (2ZZ/SZ  -  8)  +  6XZ  cos"1  <Z/S)  (14) 

la  -  B3Sa/2  -  B«/2S*  +  Y  CQ  <4  -  Tz/Sz>  -  3T  cos“MT/S>  3 

+  6XY  C  (X/T)  cos-MT/S)  -  (P/S) 2  cos~MX/P)3  (15) 

Note  that  the  two  cos_1(T/S)  terms  in  Eq.  (15)  can  be  combined; 
the  form  shown  makes  the  integration  clearer.  If  one  combines 
these  terms,  a  factor  (Xz  -  Za)  is  obtained.  As  the  limits  of 
the  third  domain  are  unchanged  by  interchanging  X  and  Z,  one  may 
sum  over  such  paired  permutations,  totaling  zero  for  these  terms. 
Therefore,  when  summed,  I3  is  equivalent  to 

I B3S3/2  -  B«/2SZ  +  YQ <4  -  Tz/Sz>  -  6XY (P/S) zcos"MX/P) .  (15E) 
for  each  permutation. 
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Fig.  1.  Chord  length  distribution  in  rectangular  volumes 
of  2^xl8xl6  (uoper)  and  20x20x10  (lower)  micrometers. 
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The  quantity  desired  in  many  calculations  is  the  integral 
■from  zero  or  to  the  maximum  length,  R,  rather  than  in  an  interior 
region.  We  shall  let  J  be  the  integral  -from  zero  to  S.  In  the 
-first  domain,  where  0  <  S  <  Z,  this  is  merely 

Jx  mJ  "n1k  dS  -  Ii<S>,  (16) 

as  Ii(O)  »  0.  In  the  second  region,  with  Z  <  S  <  T,  one  has 

Ja  -  I2(S>  -  la (Z)  +  I»(Z>  -  li(0). 


Upon  evaluation,  this  is 

Ja  -  Ia<S)  +  B*,  (17) 

where  B*.  ■  3YaZa/Ka  +  6XY  tan-MY/X).  The  integral  -from  zero  in 
the  third  domain,  where  T  <  S  <  R,  is 


j3  -  I3K (S)  +  B*  +  5Za  -  8Xa  +  6XZ  cos-1 (Z/T) .  (IB) 


It  is  easily  seen  that,  when  summed  oyer  paired  perir  ations  with 
X  and  Z  exchanged,  this  is  equivalent  to 

-  I3«(S)  +  B*  -  3Xa  +  1.5  tX.  (18E) 


-for  each  permutation. 

FKictiap  jfcaJlaiiiinufli  ..Ctiflrri-Lfinatri 

The  -fraction  o-f  chords  longer  than  S  is 


When  integrated, 


C(S)  -  1 


F  dS. 


Eq.  (1)  is  -found  to  be 


3  »  CLW  +  LH  +  WH3 (1  -  C>  -  Z  J. 


(19) 


(20) 


The  value  o-f  C  is  plotted  in  Figs.  4-6  for  the  same  dimensions  as 
employed  for  the  differential  probability  in  Figs.  1-3.  Values 
obtained  with  two  approximations  to  C,  presented  below,  are  also 
shown.  For  S  near  R,  C  is  proportional  to  (R  -  S)A. 

A  simpler  form  for  C  is  highly  useful  when  integrated  with  a 
flux  given  as  $(£),  where  £  is  the  linear  energy  transfer  (LET). 
Burke7',  Shapiro  et  al.°,  and  Petersen  et  al.*  have  reported  on 
single  event  upsets  in  fluxes  so  specified.  An  upset  will  occur 
whenever  S£  is  at  least  the  critical  energy  deposition,  E,  of 
the  device.  Therefore  we  have 

Upset  rate  ■  J*  < |»(£)  C(S)  d£,  (21) 


with  S  *  E/£.  The  upper  limit  of  integration  is  the  maximum 
value  of  £|  the  lowwr  limit  is  E/R.  If  6  and  C  are  each  given 


as  the  sum  of  one  or  more  terms  in  powers  of  £  or  S,  analytic 
integration  is  elementary. 


Approximations  £flC__£ 

Bradford10  has  developed  an  approximation  for  the  case  where 
the  minimum  dimension,  H,  is  not  more  than  a  third  of  L  or  W. 
His  formula,  utilized  in  the  above-mentioned  reports7'-’*,  is 

C  -  1  -  0.25  S/H  for  8  <  H  (22a) 

and  C  -  0.75  (H/S)a-=  for  S  >  H.  (22b) 

This  formula  is  a  fair  representation  of  the  integrated  chord 
distribution  except  for  the  longest  cases?  see  the  figures.  When 
high-LET  particles  (requiring  only  moderate  S  for  an  upset)  are 
plentiful,  the  discrepancy  at  larger  values  of  S  is  tolerable. 
When  particles  with  sufficiently  large  LET  are  few  (or  absent), 
this  approximation  cannot  be  employed. 

In  order  to  allow  use  for  spectra  where  longer  chord  lengths 
are  important,  a  similar  but  more  complex  formula  is  presented 
here.  We  first  define  a  coefficient, 

A  »  2.37/(1.80  +  R/H).  (23) 
The  two-part  formula  is 

C  -  1  -  AS/H  for  S  <  H  (24a) 

and  C  -  ( 1-A)  C  (R3- «— g»- «)  /  (R»-  *»)  ]3  (H/S) 3  for  S  >  H.  (24b) 

With  the  use  of  two  dimensions  (H  and  R> ,  the  fit  is  seen  to  be 
very  good,  being  limited  by  the  breaks  in  the  true  value.  The 
equation  has  the  pleasing  property  of  going  to  zero  at  S  ■  R. 

Let  us  consider  the  upset  rate  due  to  a  flux  proportional  to 
£~3’ ia,  with  an  abrupt  drop  to  zero  at  £M«x.  This  fits  the  high- 
LET  cosmic  ray  flux  (Figs.  1,  2,  and  4  of  Ref.  8)’  for  a  factor 
of  ten  in  LET,  but  the  abrupt  drop  is  a  factor  of  25.  When  the 
integral,  Eq.  (21),  is  transformed  into  an  integral  over  S,  the 
limits  are  SnIN  -  E/£Mfi*  to  R. 

With  this  flux,  and  for  all  six  example  geometries,  analytic 
integration  was  performed  using  the  approximation  formulas,  and 
numerical  integration  was  done  using  the  exact  equations.  Ratios 
of  "approximate  upsets"  to  the  exact  values,  as  a  function  of 
Sri  in,  are  as  expected  from  Figs.  4  to  6.  Even  for  this  extreme 
spectrum,  Bradford’s  equations  do  rather  well  in  the  region  for 
which  they  were  designed  -  flat  devices  and  small  Smin.  The  new, 
more  complex,  approximate  equation  does  better,  but  also  falters 
in  near-cubic  geometries  if  only  long  chords  can  produce  upsets. 
For  the  four  flatter  geometries,  the  upset  ratio  is  unity  within 
a  factor  of  2.42  for  Smin  <  0.5  R  using  Eqs.  (22)  and  within  a 
factor  of  1.44  for  uhin  (  0. 95  R  using  Eqs.  (24) • 


C  {PROBABILITY  OF  LENGTH  >  S) 


Fig.  Integral  chord  length  distribution 

in  the  same  volumes  as  Fig.  1. 
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C  (PROBABILITY  OF  LENGTH  >  S) 
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Fig.  5-  Integral  chord  length  distribution 
in  Che  same  volumes  as  Fig.  2. 
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Fig.  b.  Integral  chord  length  distribution 
in  the  same  volumes  as  Fig.  3- 


Summary 


The  isotropic  chord  length  distribution  is  determined  -for 
rectangular  parallelepipeds  as  a  -function  o-f  the  chord  length,  S. 
The  normalized  distribution  -function,  F,  is  given  by  Eq.  (1), 
which  uses  the  sum  over  six  permutations.  The  value  o-f  Ni«,  Na, 
and/or  N3k  is  used  as  appropriate. 

The  integral  chord  length  distribution  is  obtained  -from  the 
inde-finite  integral  <I)  o-f  N  and  the  de-finite  integral  (J)  o-f  N 
-from  zero  to  S.  The  normalized  total  integral  distribution  is 
given  in  Eq.  < 20 >. 

Results  -for  the  integrated  distribution  are  compared  with  two 
approximations,  one  o-f  them  introduced  here.  The  approximations 
are  designed  -for  integration  with  particle  -fluxes  given  in  terms 
o-f  linear  energy  transfer. 
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